Relaxor ferroelectricity induced by electron correlations in a molecular dimer Mott 

insulator 
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We have investigated the dielectric response in an antiferromagnetic dimer-Mott insula- 
tor /3'-(BEDT-TTF)2lCi2 with square lattice, compared to a spin liquid candidate k-(BEDT- 
TTF) 2 Cu2 (CN) 3 . Temperature dependence of the dielectric constant shows a peak structure obeying 
Curie- Weiss law with strong frequency dependence. We found an anisotropic ferroelectricity by py- 
rocurrent measurements, which suggests the charge disproportionation in a dimer. The ferroelectric 
actual charge freezing temperature is related to the antiferromagnetic interaction, which is expected 
to the charge-spin coupled degrees of freedom in the system. 

PACS numbers: 77.22.-d, 75.85.+t, 75.10.Kt 
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The discovery of new types of multiferroics [l], 
where ferroelectricity is driven by nontrivial magnetism, 
have also initiated new research into a more fundamen- 
tal understanding of spin-charge coupled phenomena in 
strongly correlated systems. However, controlling fer- 
roelectricity is difficult in such systems due to strong 
spin structure-electrical dipole coupling. Thus, electronic 
ferroelectricity 0, 0] is now explored as an alternative 
route towards multiferroics, related to the charge or- 
dered state in Mott insulators 0, Among them, 
molecular dimer Mott insulators are a good candidate 
for the ferroelectricity, which are typical strongly corre- 
lated electron systems, constructed by dimers of organic 
molecules such as bis(ethylenedithio)tetrathiafulvalene 
(BEDT-TTF). When each dimer in a crystal possesses 
a charge carrier, in particular for the half- filled case in a 
dimer molecular orbital, the electron system often shows 
a Mott transition for relatively large (on-site) Coulomb 
interactions, U , at the dimer site compared to the elec- 
tron hopping energy, t Q . An important physical insight 
is that a dimer, rather than a molecule, can be viewed 
as the primitive constituent in a crystal. In this way, 
the understanding of strongly correlated electron systems 
based on the dimer picture within the Hubbard model 
has progressed 0, H|. However, recent reports on dielec- 
tric properties [9( imply that the charge degrees of free- 
dom in a dimer survive even in the dimer Mott insulating 
phase T3, 11 1 . This discovery will enable a deeper under- 
standing of organic Mott insulators beyond the successful 
dimer lattice model and the discovery of new function- 
ality. The origin of this ferroelectricity is different from 
the recently discovered hydrogen-bonded organic ferro- 
electrics [l2j and the Peierls instability-induced inversion 
symmetry breaking in one-dimensional systems [HI . 

An anomaly in the dielectric constant has been 
observed for a dimer Mott insulator «-(BEDT- 
TTF) 2 Cu 2 (CN) 3 i with a triangular lattice in terms 



of the dimer picture, which is one candidate for a spin 
liquid 14-1^. The temperature dependence of the di- 
electric constant for k-(BEDT-TTF) 2 Cu 2 (CN) 3 shows 
a ferroelectric relaxor-like frequency dependence, ex- 
pected to exhibit charge disproportionation [ll| in a 
dimer or st rong charge fluctuations near a charge or- 
dered state [l0|. The dielectric constant in k-(BEDT- 
TTF) 2 Cu 2 (CN) 3 obeys the Curie Weiss law with a Curie 
temperature of 6 K. However, there is no explicit evi- 
dence of ferroelectricity or the charge freezing typically 
seen alongside relaxor behaviour. This dielectric anomaly 
may be related to the quantum spin liquid nature of this 
material, which further complicates the physical origin 
of the dielectric constant anomaly. Here, we have ob- 
served a similar dielectric anomaly due to the response of 
charge degrees of freedom in a dimer for a simpler dimer 
Mott insulator /3'-(BEDT-TTF) 2 ICl 2 with a square lat- 
tice and discovered evidence for relaxor ferroelectricity by 
pyrocurrent measurements with an anisotropy consistent 
with the dimer arrangement. These results may lead to 
the future discovery of electronic ferroelectric materials 
in dimer Mott insulators and a deeper understanding of 
the spin liquid state. 

/3'-(BEDT-TTF) 2 ICl 2 is an anisotropic quasi-two di- 
mensional molecular dimer Mott insulator with one hole 
carrier in the dimer of a BEDT-TTF molecule. The unit 
cell of /3'-(BEDT-TTF) 2 ICl 2 is a square lattice, in con- 
trast to k-(BEDT-TTF) 2 Cu 2 (CN) 3 . The space group of 
this system is triclinic PI where the inversion symme- 
try produces one dimer in a unit cell, as shown in Fig . 
1(a) [13] • Dimerization was also supported by NMR flp| 
and a calculation of the intradimer transfer integral jl7j | . 
Thus, the dimer picture is applicable for this system, 
and this material is considered a half-filled Mott insula- 
tor with an antiferromagnetic transition at Tn = 22 K 
[l9j | . The higher magnetic transition temperature indi- 
cates that this system is located within the well-defined 
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FIG. 1: (Color online) Crystal structure of /3'-(BEDT- 
TTF) 2 IC1 2 . (a) BEDT-TTF molecule arrangement viewed 
along the long axis of BEDT-TTF molecules and (b) the crys- 
tal structure seen from the c axis. The blue ovals represent 
dimerized molecules symmetric with respect to the inversion 
operation at the point represented by open circles. The unit 
cell is represented by blue lines in (b). Hydrogen atoms are 
not drawn for clarity. 



dimer Mott insulating phase [6j, in comparison with k- 
(BEDT-TTF) 2 Cu 2 (CN) 3 , which docs not exhibit long 
range magnetic order down to at least 32 mK 14|. As 
shown in Fig. 1. one dimer exists in a unit cell, aligned 
along the b axis 1JJ . The dimers are arrayed in a column 
on the be plane, and the layers are stacked along the a* 
axis (_L be plane), well separated by ICI2 anion layers. 
Band structure calculations [l7| and optical conductivity 
measurements [20I , [2l| indicated strong one-dimensional 
anisotropy along the b axis (in the be plane). The sys- 
tem becomes superconducting under 8 GPapressure 17| , 
as typically seen in dimer Mott systems [8j. Thus, this 
system could be understood as a typical dimcr-Mott insu- 
lator without considering the presence of ferroelectricity. 

The samples used in this study were grown by an elec- 
trochemical method [22J • Typical dimensions of the spec- 
imens are 0.7 mm x 0.7 mm x 2 mm along the a* , b* , and 
c axes, respectively. An electric field of 5 V/cm was ap- 
plied along different crystallographic directions in order 
to measure the dielectric constant with an LCR meter. 
Pyrocurrent measurements were performed with increas- 
ing temperature at a rate of 10 K/min after applying a 
poling electric field from -1.2 kV/cm to 4.8 kV/cm at 
approximately 100 K. 

Figure 2 shows the temperature dependence of the di- 
electric constant e for /3'-(BEDT-TTF) 2 ICl 2 in (a) E || a* 
perpendicular to the two dimensional BEDT-TTF layer 
and (b) E \\ b parallel to the dimer array at each fre- 
quency (20 to 200 kHz). An anomaly with a strong 
frequency dependence is discerned in both directions at 
around 80-150 K. The anomalies at high frequencies, e.g. 
> 50 kHz, show up at 130-150 K. With decreasing fre- 
quency, the anomaly in the dielectric constant becomes 
large and forms a peak structure. A similar frequency 
dependent dielectric constant has been observed for k- 
(BEDT-TTF) 2 Cu 2 (CN) 3 , which is typically seen in rc- 
laxor ferroelectric materials with randomness and frus- 




50 100 
Temperature (K) 



FIG. 2: (Color online) Temperature dependence of the di- 
electric constant, e, in (a) E \\ a* (± be plane) at frequencies 
from 500 Hz to 100 kHz and (b) E \\ b from 20 Hz to 200 
kHz for /3'-(BEDT-TTF) 2 ICl 2 . The data in the high tem- 
perature region was omitted for clarity. The inset shows the 
Vogel-Fulcher (VF) relation using the dielectric constant at 
frequencies from 20 Hz to 7 kHz indicated as open triangles 
in (b). The black solid line in the inset is a best fit to the VF 
relation. T max and Typ = 23 K represent the temperatures 
at which e takes a maximum value at each frequency and the 
charge freezing temperature obtained by the VF fitting, re- 
spectively. Solid triangles in (b) represent the expected T max 
by the relation indicated the filled circles in the inset. 



tration. In /3'-(BEDT-TTF) 2 ICl 2 at 500 Hz, the peak 
value from the background level at the lowest temper- 
ature corresponds to a Ae = 10 which is also compa- 
rable to that in k-(BEDT-TTF) 2 Cu 2 (CN) 3 Q at the 
corresponding frequency and crystallographic direction. 
While we could not obtain reliable data for e at frequen- 
cies lower than 500 Hz in E \\ a* , we were able to measure 
E || b. At the lowest frequency of 20 Hz, the anomaly 
grows as large as Ae = 20 for E || b. Thus, dielec- 
tric anomalies discerned in the both directions arc of an 
identical origin, which is expected as charge dispropor- 
tionation occurs in a dimer along the b axis. 

The inverse dielectric constant subtracted by a con- 
stant is shown in Fig. 3(a), which clearly obeys the 
Curie- Weiss law with a constant term e CO nst.: e/^o = 
T _ Tc + Econst./^o- Here eo, C, and Tq are the vacuum 
permittivity, Curie constant, and Curie temperature, re- 
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spectively. The black line in Fig. 3(a) represents the 
best fit obtained for the data at frequencies less than 
5 kHz. The constant term e CO nst./eo is 13.8. Accord- 
ing to the slope of the inverse dielectric constant, C, 
the dipole moment parallel to the b axis is estimated 
to be 0.13ed, where e and d are the electron charge 
and distance between the BEDT-TTF molecules in a 
dimer, 3.6 respectively. The Curie temperature 

Tq was found to be 67 K. This Tq is ten times larger 
than that for ^-(BEDT-TTF^CWCN^, which is un- 
derstood in terms of a large U/t |6j. Since /?'- (BEDT- 
TTF) 2 ICI2 is located within the well-defined antiferro- 
magnetic Mott insulating region, the energy scale for an 
electronic phase transition should be higher than that 
of k-(BEDT-TTF) 2 Cu 2 (CN) 3 , especially considering the 
effects of geometrical frustration. 

While the apparent ferroelectric Curie temperature Xc 
and relaxor-like charge-freezing or Vogel-Fulcher tem- 
perature, Xvf, of the domains reported for k-(BEDT- 
TTF) 2 Cu 2 (CN) 3 are identical at 6 K 0, these tempera- 
tures for /3'-(BEDT-TTF) 2 ICl 2 , as described below, are 
clearly different as shown in Fig. 3(a) and the inset of Fig. 
2, respectively. The inset of Fig. 2 shows a fit to the data 
for E || b of /3'-(BEDT-TTF) 2 ICl 2 to the Vogel-Fulcher 
(VF) relation for e: f = f exp {-E VF /(T max - T VF )} 
where /o, Tvf, and T max , are the representative flipping 
frequency of electric dipoles, average random energy bar- 
rier, and temperature at which e takes a maximum, re- 
spectively While fitting to this relation, e as measured 
at frequencies from 20 Hz to 7 kHz in Fig. 2(b) was 
used. The empirical VF relation well represents the fre- 
quency dependence of the peak in e with a Typ = 23 
K, E VF = 804 K, and f = 2.17 x f0 7 Hz. In this sys- 
tem, the charge freezing temperature, TVf, is much dif- 
ferent from the Curie temperature Tq = 67 K, although 
these temperatures arc often considered to be the same. 
Since the VF relation reveals little about the temperature 
dependence of the dielectric constant, we measured the 
pyrocurrent to uncover the actual charge freezing tem- 
perature. 

Figure 3(b) shows the electric polarization, P, with 
E || b obtained by pyrocurrent measurements and (c) 
the temperature dependence of the magnetic suscepti- 
bility, Xi i R a magnetic field with H || a*, b, and c for 
/3'-(BEDT-TTF) 2 ICl 2 . A small but measurable polar- 
ization was observed with a critical temperature of Tfe 
= 62 K, which is nearly identical to the estimated Tq = 
67 K. P increases linearly with poling electric field, E po i, 
at least up to 2 kV/cm. Leakage currents prevented the 
measurement of P for |£p i| larger than 5 kV/cm. The 
electric field linearity and small polarization values in- 
dicate that the fully polarized state has not yet been 
obtained under poling electric fields of 5 kV/cm. The 
expected full moment is estimated as 2 //C/cm 2 assum- 
ing the dipole moment of 0.13ed, which is comparable to 
hydrogen bonded ferroelectric materials (l2j . However, 
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FIG. 3: (Color online) (a) Temperature dependence of in- 
verse dielectric constant at each frequency in /3'-(BEDT- 
TTF)2lCl2- A constant term of e C onst./eo — 13.8 was sub- 
tracted from the data. The black straight line represents the 
best fit result to the Curie- Weiss law. Temperature depen- 
dence of (b) polarization with E \\ b and a poling electric 
field from -1.211 to +4.844 kV/cm obtained by pyrocurrent 
measurements, and (c) magnetic susceptibility with H \\ a* , 
b and c for /3'-(BEDT-TTF) 2 ICl 2 . The polarization at 7 K 
for each poling field in (b) is plotted in the inset, which does 
not correspond to the standard P — E curve at 7 K. The 
olarization exhibits a tiny anomaly at Tn. The susceptibility 
clearly shows an anisotropy at the antiferromagnetic transi- 
tion temperature, Tn, but no anomaly near the ferroelectric 
temperature, Tfe. 



the experimentally obtained polarization of 13 nC/cm 2 
at .Epoi = 4.8 kV/cm is small. Thus, such a polarization 
indicates the formation of glassy polar domains, also con- 
sistent with the strong frequency dependence in e (Figs. 
2(a) and 2(b)). A pyrocurrent while poling with E || a* 
could not be observed up to 3 kV/cm, which supports the 
assertion that electric dipoles align along the dimer ar- 
rangement. The temperature dependence of \ (see Fig. 
3(c)) does not exhibit a notable anomaly at Tfe, per- 
haps because the glassy ferroelectricity observed here is 
only electric field-induced. However, a small anomaly in 



the P-T curve is seen around Tn. The magnetic phase 
transition at Tn is thermodynamically well defined and 
may affect the electric dipoles different from the glassy 
transition. Therefore, in this system, the spin and charge 
degrees of freedom within a dimer have a moderate cou- 
pling, even though the antiferromagnetic and ferroelec- 
tric transition temperatures are different. 

Our discovery of ferroelectricity with a strong fre- 
quency dependence, similar to that in «-(BEDT- 
TTF) 2 Cu2(CN) 3 , may be due to the intrinsic phe- 
nomenon of electronic ferroelectricity in a dimer-Mott 
insulator with antiferromagnetic interactions. Although 
the mechanism for slow dynamics with possible domain 
formation is an open issue, the discrepancy between the 
actual charge freezing Tfe (— Tc) and the fit result of 
Tvf originates from spin-charge coupling as shown be- 
low. The charge distribution in the dimer-Mott insulat- 
ing state has the centre of inversion as shown in Fig. 4(a). 
However, ferroelectricity or charge disproportionation in 
a dimer is driven intrinsically by the hybridization of two 
orbitals with different parity [231 ] . Such partial hybridiza- 
tion produces electrically polar excited states, considered 
to result from collective phenomena induced by the inter 
dimer Coulomb correlation, V (Fig. 4(b) and 4(c)). This 
interaction can lead to char ge fl uctuations in a dimer, as 
predicted theoretically lllf . showing a Curie- Weiss 
behaviour for the dielectric constant. However, when 
the temperature of the system decreases below the en- 
er gy s cale of the antiferromagnetic interaction Jaf (= 59 
K [l9(), the approach of electrons with opposite spins is 
favourable due to the increase in transfer as Jaf = —t 2 /U 
(Fig . 4(d)), which may be relevant to the increase in Tn 
24 1 and the presence of superconductivity under pressure 
17 1. Thus, an antifcrroelectric interaction is possibly re- 
alized with Jaf rather than a ferroelectric interaction 
without Jaf- In this sense, Tvf is almost the same as Tn 
and is related to the antiferroelectric interaction between 
ferroelectric domains. Therefore, additional randomness 
and frustration between ferroelectric and antiferroelec- 
tric interactions are produced by the spin-charge coupled 
degrees of freedom, resulting in a charge freezing at Tc 
higher than Tvf- 

Figure 4(e) shows a schematic picture of the electronic 
phases with respect to the parameters Tvf and Tq in 
/?'-(BEDT-TTF) 2 ICl 2 and k-(BEDT-TTF) 2 Cu 2 (CN) 3 . 
The discrepancy in the phases and parameters is clearly 
seen between the crystals. In the case of k-(BEDT- 
TTF) 2 Cu 2 (CN) 3 , T c and T VF arc the same, imply- 
ing that the competition between ferroelectric and 
antiferroelectric interactions, although magnetic and 
electric phase transitions have not been observed at 
any temperature. This is contrastive to «-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Cl [H] which was reported to show 
a possible ferroelectricity. The frustration amongst 
spins and charges between dimers may be responsi- 
ble for the quantum spin liquid nature of k-(BEDT- 
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FIG. 4: (Color online) Schematic pictures of (a) dimer-Mott 
insulator, inter dimer Coulomb correlation (V) with (b) ferro- 
electric and (c) antiferroelectric dimers, and (d) antiferromag- 
netic interaction (Jaf) with antiferroelectric dimers within 
the be layer, where a dimer is represented by an open circle. 
Solid circles and arrows represent charge and spin, respec- 
tively, (e) Schematic electronic phase diagram for /3'-(BEDT- 
TTF) 2 IC1 2 and k-(BEDT-TTF) 2 Cu 2 (CN) 3 , where the dielec- 
tric constant starts to increase at around 120 K and 60 K, 
respectively. Tn and Tfe correspond to antiferromagnetic 
and ferroelectric transition temperatures, respectively. The 
change in parameters of the Curie- Weiss (Tc) and the Vogel- 
Fulcher (Tvf) temperatures are also shown. 



TTF) 2 Cu 2 (CN) 3 with a unique structure. 

In conclusion, we have revealed the dielectric re- 
sponse in the antiferromagnetic dimer-Mott insulator 
/3'-(BEDT-TTF) 2 iCl 2 with a square lattice, compared 
to the spin liquid candidate k-(BEDT-TTF) 2 Cu 2 (CN) 3 . 
We observed a peak structure in e with strong frequency 
dependence, probably due to ferroelectric relaxor-likc 
domain formation. The anisotropic pyrocurrent corre- 
sponds closely with charge disproportionation in a dimer, 
and shows the actual charge freezing temperature is re- 
lated to the antiferromagnetic interaction. This, in turn, 
supports the finite spin-charge coupled degrees of free- 
dom in the system. Our result confirms the importance of 
the charge degrees of freedom, moderately coupled with 
spin, in BEDT-TTF dimers, as a key concept widely 
applicable to molecular dimer Mott transition systems 
and toward the discovery of electronic ferroelectricity in 
dimer Mott insulators. 

The authors would like to thank S. Ishihara, S. Iwai, 
and I. Terasaki for valuable discussions, and acknowl- 
edge the experimental contributions by M. Abdcl-Jawad 



5 



at the initial stage of this study. This work was sup- 
ported by Grants-in-Aid for Scientific Research from 
JSPS (24540357) and MEXT, Japan (23110702 and 
23102502). 



[1] T. Kimura, T. Goto, H. Shintani, H. Ishizaka, T. Arima, 

and Y. Tokura, Nature 426, 55 (2003). 
[2] H. Katsura, N. Nagaosa, and A. V. Balatsky, Phys. Rev. 

Lett. 95, 057205 (2005). 
[3] J. van den Brink, and D. Khomskii, J. Phys.: Condens. 

Matter 20, 434217 (2008). 
[4] S. Ishihara, J. Phys. Soc. Jpn. 79, 011010 (2010). 
[5] M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 

70, 1039 (1998). 
[6] K. Kanoda, J. Phys. Soc. Jpn. 75, 051007 (2006). 
[7] H. Kino, and H. Fukuyama, J. Phys. Soc. Jpn. 64, 2726 

(1995). 

[8] H. Seo, C. Hotta, and H. Fukuyama, Chem. Rev. 104, 
5005 (2004). 

[9] M. Abdel-Jawad, I. Terasaki, T. Sasaki, N. Yoneyama, 
N. Kobayashi, Y. Uesu, and C. Hotta, Phys. Rev. B 82, 
125119 (2010). 

[10] M. Naka, and S. Ishihara, J. Phys. Soc. Jpn. 79, 063707 
(2010). 

[11] C. Hotta, Phys. Rev. B 82, 241104(R) (2010). 
[12] S. Horiuchi, and Y. Tokura, Nature Mater. 7, 357 (2008). 
[13] F. Nad, P. Monceau, C. Carcel, and J. M. Fabre, Phys. 
Rev. B 62, 1753 (2000). 



[14] Y. Shimizu, K. Miyagawa, K. Kanoda, M. Maesato, and 

G. Saito, Phys. Rev. Lett. 91, 107001 (2003). 
[15] S. Yamashita, Y. Nakazawa, M. Oguni, Y. Oshima, H. 

Nojiri, Y. Shimizu, K. Miyagawa, and K. Kanoda, Nature 

Phys. 4, 459 (2008). 
[16] M. Yamashita, N. Nakata, Y. Kasahara, T. Sasaki, N. 

Yoneyama, N. Kobayashi, S. Fujimoto, T. Shibauchi, and 

Y. Matsuda, Nature Phys. 5, 44 (2008). 
[17] H. Kobayashi, R. Kato, A. Kobayashi, G. Saito, M. Toku- 

moto, H. Anzai, and T. Ishiguro, Chem. Lett. 15, 89 

(1986). 

[18] K. Miyagawa, K. Kanoda, and A. Kawamoto, Chem. 

Rev. 104, 5635 (2004). 
[19] N. Yoneyama, A. Miyazaki, T. Enoki, and G. Saito, Bull. 

Chem. Soc. Jpn. 72, 639 (1999). 
[20] H. Kuroda, K. Yakushi, H. Tajima, A. Ugawa, M. 

Tamura, Y. Okawa, A. Kobayashi, R. Kato, H. 

Kobayashi, and G. Saito, Synth. Met. 27, A491 (1988). 
[21] H. Taniguchi, M. Miyashita, K. Uchiyama, K. Satoh, N. 

Mori, H. Okamoto, K. Miyagawa, K. Kanoda, M. Hedo, 

and Y. Uwatoko, J. Phys. Soc. Jpn. 72, 468 (2003). 
[22] H. Taniguchi, R. Sato, K. Satoh, A. Kawamoto, H. 

Okamoto, T. Kobayasi, and K. Mizuno, J. Low Temp. 

Phys. 142, 437 (2006). 
[23] S. Onoda, S. Murakami, and N. Nagaosa, Phys. Rev. 

Lett. 93, 167602 (2004). 
[24] Y. Eto, and A. Kawamoto, Phys. Rev. B 81, 020512(R) 

(2010). 

[25] P. Lunkenheimer, J. Miiller, S. Krohns, F. Schrettle, A. 
Loidl, B. Hartmann, R. Rommel, M. de Souza, C. Hotta, 
J.A. Schlueter, and M. Lang. larXiv:1111.2752V 2 (2011). 



